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Abstract 



A model of unification of dark matter and dark energy based on the modeling 
of the speed of sound as a function of the parameter of the equation of state 
is introduced. It is found that the model in which the speed of sound depends 
on the power of the parameter of the equation of state, c 2 s = a{— w) 1 , contains 
the generalized Chaplygin gas models as its subclass. An effective scalar field 
description of the model is obtained in a parametric form which in some cases 
can be translated into a closed form solution for the scalar field potential. A 
constraint on model parameters is obtained using the observational data on the 
Hubble parameter at different redshifts. 



1 Introduction 

The expansion of the universe is one of the most fascinating phenomena that sci- 
ence has encountered so far. It has served as a rich source of information on the 
nature and the composition of the universe. One of recently established aston- 
ishing features of cosmic expansion is that it is currently undergoing a phase of 
acceleration [TJ |2j [3j 0] . The source of this acceleration has not yet been unam- 
biguously identified, although many proposals for its nature have been put forward 
(see [6] and references therein) . The existence of dark energy, a mysterious com- 
ponent with negative pressure is still the most serious candidate. Another dark 
component of the universe, dark matter, seems to leave its imprint on astrophysical 
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and cosmological scales, ranging from galaxies to galactic clusters and large scale 
structure in the universe. 

The idea that both dark matter and dark energy are actually the manifes- 
tations of a single dark component is both natural and appealing. It appeared 
early in the literature and its the most acclaimed representative is probably the 
Chaplygin gas [7] as a model of unification of dark matter and dark energy [8]. 
The class of unifying DM-DE models is often referred to as quartessence QH] . 
This phenomenologically introduced model can be motivated from string theory 
[E1QT1III]. Its usually studied extension, the generalized Chaplygin gas model, was 
first introduced in [13]. The agreement of the Chaplygin gas and its extensions 
with observations has been extensively tested, including the analyses with the su- 
pernova la data [21 H5] . CMB [IB] , observable Hubble parameter data [TTJ, [18] 
and large scale structure observations [191 EQl [ZD 1221 [23] , including the nonlinear 
evolution in structure formation [21] and gravastar formation [2S]. Different data 
can be combined to produce tighter parameter constrains such as in |26U271 l28 l l29|. 
Strong constraints on the generalized Chaplygin gas have been obtained that ques- 
tion its viability as a cosmological model distinguishable from the ACDM model. 
In order to better accommodate observational constrains, various unified models 
based on Chaplygin gas have been proposed, such as the modified Chaplygin gas 
model [301 EI] ■> recently reviewed and constrained in [32], [331 El] or hybrid Chaply- 
gin gas leading to transient acceleration [35]. The fact that perfect fluid model can 
be fully described by defining speed of sound equation has been used in [36] . The 
idea of DM-DE unification with non-canonical scalar fields has been recently stud- 
ied in [37, 38j. An interesting model called Dusty Dark Energy, recently introduced 
in [35], achieves the DM-DE unification in the formalism of the X(p-fluid, resulting 
in the zero speed of sound and one scalar degree of freedom. Other approaches to 
models of DM-DE unification that avoid the speed of sound problem are purely 
kinetic k-essence models [4Q |. |4"T1 |4"2] and tachyon models [431 144] . 

The structure of the paper is the following. After the introduction presented in 
this section, in the second section a general class of barotropic fluid models defined 
by the function c 2 . is discussed. In the third section the model with the constant 
speed of sound is studied and in the fourth section the principal model of the 
paper, defined by c 2 s {w) = a(— u>) 7 is introduced. The fifth section is focused on an 
effective representation of the model in terms of a minimally coupled scalar field. In 
the sixth section the comparison of the model prediction against the observational 
data on the Hubble parameter at different redshifts is made and the constraints 
on the model parameters are presented. The seventh section closes the paper with 
the discussion and conclusions. The Appendix outlines an approach in which the 
solution with the piecewise constant speed of sound is used as an approximation 
of the dynamics of a fluid with a general dependence of c 2 s on w. 



2 The model 

The equation of state of a barotropic cosmic fluid can in general be written as an 
implicitly defined relation between the fluid pressure p and its energy density p 

F(p,p) = 0. (1) 
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The parameter of the equation of state w = p/p can be used to substitute the 
pressure, p = wp, so that the relation ([I]) becomes G(p, w) = 0, with G(p, w) = 
F(p,p). This relation implies that p and p can be considered as functions of w, i.e. 
p = p{w) and p = p(w) = wp(w). 

The speed of sound waves in the barotropic cosmic fluid is defined as 



From (H|) it follows 



which leads to 



9F , dF , 

- dp+ - dp = 0, (3) 



dF 
dp 



Inserting the relation p = wp into (|3|) and using @, one readily obtains 

dp dw 



p (? a — w 



(5) 



Combining this expression with the continuity equation for the fluid 



one finally obtains 



dp + 3p(l + w)— =0, (6) 
a 



dw da dz 

-3— = 3- . (7) 



(cf — w)(l + w) a 1 + z 

As p and p are functions of w, the expression for the speed of sound can be written 
as 

c 2 = ± = jt (8) 

s dp * ' w 

i.e. = Cg(io). Therefore, the expression ([7|) is a dynamical law for the EOS 
parameter w. 

The only unknown part in ([7]) is the functional form of c^(w). By its very 
definition there are some constraints on it, such as that it should be nonnegative 
and smaller than the speed of light squared, c 2 . The modeling of cosmic fluid 
unifying dark matter and dark energy by modeling (? s as a function of w may, 
therefore, be more suitable since the requirements on c 2 can be immediately built 
in. Furthermore, the properties of the solutions of ([7]) depend on the zeros of the 
function <? s (w) — w. In particular, the dynamics of w(z) is confined to intervals 
determined by the zeros of c^(w) — w and w = —1. For some models with a variable 
speed of sound see [351 H6| 137 } HH] . 

In the following two sections we discuss specific choices for (? s and discuss their 
relation with models previously studied in the literature. 
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3 Constant speed of sound 



The simplest possibility is a constant speed of sound. This possibility was recently 
studied in , see also [H2 [SDJ EH E2] . A direct inspection of © shows that for 
the case of constant speed of sound w is confined to one of the intervals: (—00, —1), 



(-1, C 2) and (cj 
factor as 



00 . 



For c 2 = const, the parameter of EOS evolves with the scale 

|3(l+cJ) 



w 



C 2 4+ML(1 + Z ) 



1 



1 + WQ 



(1 + Z )3(l+c2) + 1 



with w(0) 



and 



wq. From this relation it immediately follows 

w 1 + Wq 



P = Po- 



ll 



C 2 sP 



Po{c 2 s 



Wq 
- W 



w ) 



Po- 



Po- 



1 + C 2 



w 



-(1 + z 



3(1+4) 



+ 1 



c 2l±^L (1 + z) 3(l + ^) 



c% - Wo 



Finally we obtain 



P + P = Po(l + wo){l + z 



,3(l+c2) 



(9) 



(10) 



(11) 



(12) 



For — 1 < w < Cg, the energy density remains positive throughout the evolution 
of the universe. With the cosmic expansion, the pressure evolves from positive to 
negative values. In particular, the pressure changes sign at 1 + z = {c 2 s {l+wo) / (c 2 s — 
u;o)) _1 / ( ' 3( ' 1+Cs ''* ) • For wq < — 1 the pressure is negative throughout the evolution of 
the universe, p < 0. The energy density evolves from negative to positive values 
with the expansion and it changes sign at 1 + z = (— (l + wo)/(c 2 — wq))^ 1 ^ 3 ^ 12 ^ . 
For wo > c 2 , the pressure is always positive, whereas the energy density evolves 
from positive to negative values with the expansion, with a change of sign at 

1 + Z = ((1 + wo) /{Wo - C 2))-V(3(l+a). 

The solution for c 2 = const can be used as a building block for approximating 



general functional forms c s 



w 



The discussion of a piecewise constant approxima- 



tion of the general ci redshift dependent function is given in the Appendix. 



4 Power law dependence: c 2 s = a(-w) 1 

For a more general parametrization 

c 2 = a{-wy , (13) 

the equation ([7]) in general needs to be solved numerically to obtain w = w(z). An 
example of such a numerical solution is presented in Fig[TJ The redshift dependence 
of the speed of sound is presented in Fig. [2j We are primarily interested in the 
solutions in which w is confined to the interval (—1,0) and the model is properly 
defined for w < 0. For 7 > 1 and a > the dynamics of w is confined to one of the 
intervals (— 00,— 1), (—1,0), (0,oo). For 7 > 1 and a < 0, Eq. (|7|) can be written 
as 

(1 + Z )^L = -3aw(w + m-wy- 1 - (-w.y- 1 ) , (14) 
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where to* = — (— l/a) 1 ^ 7-1 ) is the additional zero of the denominator in ([7]). For 
a < — 1, falls in the interval (—1,0) which is of primary interest in this paper. 
Depending on the relation of wq and w*, there are two distinct cases. For wo > w#, 
with the expansion w(z) evolves from w = at early times towards w* which 
it asymptotically reaches in the infinity. For wq < w* the expansion starts from 
w = — 1 and asymptotically reaches w* in the infinity. In the latter case the model 
cannot serve as a model of unification of dark matter and dark energy, but it could 
serve as a model of dark energy only. 
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Figure 1: The redshift dependence of the EOS parameter of the unified DM-DE compo- 
nent for several values of (a, 7) parameters. For all parameter values the present value 
of the EOS parameter, wq, is selected to match the value of the spatially flat ACDM 
model. 

The EOS of the cosmic fluid with this parametrization of the speed of sound, 
however, can be calculated analytically. From the definition 



dp f P 
dp V P 



(15) 



a + (-Wq) 



1-7 



a + (-wy-~f _ 



1/(1-7) 



(16) 



for 7/1, one readily obtains 

P = Po 

For 7=1, the EOS becomes 

p = Ap~ a . (17) 

The result (|17p reveals that for 7 = 1, the parametrization (|13|) is equivalent 
to generalized Chaplygin gas, whereas for 7 = 1, a = 1 the EOS of Chaplygin gas 



. 8 
. 6 



^0.4 



i i 

\ ] 

i 
i 
i 
i 
i 

_ i 

: i 
• i 



i t 



0.2 
. 







1 





ff = 2,y = 1 






a = 2,y = 3 






a = 2,y =5 






ACDM 





3 



z 



Figure 2: The redshift dependence of the speed of sound squared, c^, of the unified 
DM-DE component for several values of (a, 7) parameters. For all parameter values the 
present value of the EOS parameter, w , is selected to match the value of the spatially 
flat ACDM model. 
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is reproduced. Therefore, (|13j) represents a broad new class of unification models 
which contains the Chaplygin gas as its special case and the generalized Chaplygin 
gas as its subclass. The study of the model (fT3J) could potentially lead to deeper 
understanding of strong constraints appearing in the comparison of the Chaplygin 
gas and the generalized Chaplygin gas with the observational data. 

For 7 = 2 and a / -1, Eq. ([7]) can be solved analytically and presented in a 
closed form: 

/ i 1 \ l/(«+l) / 1 / \ a/(a+l) 

Furthermore, for a special case a = 1, explicit expressions for the evolution can 
be obtained. The function w(z) then yields 

w(z) = 1 , (19) 

l + i^(l + *)6 

whereas the energy density and pressure read 



1 \ ' ^c>-*) ( 



and 



p(z) = -wqPq — (20) 

-w + 1 



1 + 



p(z) = p wo — , (21) 

-wq + 1 

respectively. 

Finally, for a = —1/2 Eq. ()18|) can also be presented as an explicit expression 
for w{z): 

w(z) = -l + -= 1 =■ (22) 

The energy density and pressure are given by the expressions 



and 



= »77^ I 1 + , . , = I ■ (24) 



(W0 + 1) 2 

5 Scalar field representation 

If the cosmic fluid with the speed of sound (|13p is the dominant component in the 
universe, in a spatially flat FLRW universe the Hubble parameter is 

n 8vrG . . 

H 2 = ^—p. (25) 
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In this section we focus on finding the effective description of the cosmic dynamics 
with the cosmic fluid defined by (|13p in terms of a minimally coupled scalar field ip 
with a potential V(ip). For an effective representation of the cosmic fluid in terms 
of the scalar field, the expressions for p and p are 

p = \? + V(tp) (26) 

and 

p=l<j?-V(tp). (27) 

The time derivative of the scalar field can be expressed as ip 2 = p + p = 
(1 + w)p(w) from ([26]) and ([27]) . On the other hand 



Combining these two expressions for <p and using (|25p . the following equation is 
obtained 

8vrG ,11, 
3 o vl + u;(c~ — io ) 

Integration of this equation leads to (p = (p(w). In a similar way, combining (|26p 



and (|27p the scalar field potential is 



V{<p) = ~(l-w)p(w). (30) 



Therefore, assuming that ()29p is analytically integrable, the scalar field represen- 
tation is available in the parametric form as <p = <p(w) and V(<p) = V(w). Even 
in the case when it is not possible to analytically integrate (|29p . Eqs (j29[) and (j30|) 
provide a direct procedure (including numerical integration for each value of w) 
for the parametrically defined pair of quantities (p(w),V(w). In the remainder of 
this section we consider an analytical scalar field reconstruction for several specific 

examples and use the abbreviation <p = J ^py?- 
5.1 cj, = a 

For the case of constant sound speed with a > — 1, the equations (fTU|) and ([30]) 
yield 

V(<P) = \p (a-w )^^. (31) 
2 a — w 



On the other hand, Eq. (|29j) is readily integrated to 

0-0o = ^^ fin ^ + V^ _ ln v^-VTT^ (32) 
The parameter of EOS is then 

™ = (l + «)(j^) 2 -l, (33) 
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where 



y 



Vi + w + VTT 



« e =F3VI+S(^-^o) 



3 =F3a/T+^(0-0) 



VI + w - \/T-Fa 
Combining (f3Tj) . (f33|) and (f34"j) the scalar field potential is 



^0) = \po{a - w ) 



1 + 



— (cosh ( - 

l + a\ y. 



VTT 



o 



(34) 



(35) 



5.2 ci 



-aw 



In this section we consider the sound speed linearly dependent on the EOS pa- 
rameter. As already shown in section [H such a dependence is characteristic of the 
generalized Chaplygin gas. In this subsection we reconstruct the scalar potential 
for the generalized Chaplygin gas using the method of parametric definition of 4> 
and V{4>) in terms of w. The integration of (|29p gives 



±3(l + a)(0 
which results in 



In 



In 



VI + w + 1 



w 



VI + w - 1 VI + w + 1 ' 



l+y 



where 

VI + wp - 1 
VI + + l' 
The scalar field potential is 



y 



( ±3(l+a)(0-<fc>) 



3 ±3(l+a)(0- 



^) = ip (l-^)(- N ) 
2 V^o/ 



-l/(l+a) 



(36) 

(37) 
(38) 

(39) 



The combination of (|3"7|) . ([55]) and (|3"9"1) results in the potential 



2 \ w 



-l/(l+o) 



r. 



cosh I ^(1 + a)(<t>- (f>) 



-2a/(l+a) 



2/(l+a) 



cosh ( -(1 + q)(0 



(40) 



for the generalized Chaplygin model. Specifying a = 1 we obtain the potential for 
the Chaplygin gas 



1 



V{(t>) = ^(-wo) 1 ' 2 cosh(3(0 - $)) + 



1 



cosh(3(0 — cj))) 



(41) 
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5.3 



For 7 = 2 the integration of 



leads to 



V^o + 1 + 1 \ZuT+I - 1 V + P J VV^o + l-/3 



where /3 = y^(a + 1)/q. The scalar field potential has the form 

(1 — w)(aw — 1) 



V(<P) = TjPO 

2 awQ - 



1 



1//3- 



(42) 



(43) 



Equations (|42|) and (|43|) provide the solution for the scalar field potential in a 
parametric form. For 7 > 1 the evolution of the EOS parameter w is confined to 
the interval (—1, 0). Taking the values for w from this interval the plots of V{4>) are 
readily obtained. The functional forms of V(4>) for several values of a are depicted 
in Fig. El 
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Figure 3: The effective scalar field potential for several values of parameters (a, 7). The 
potential V(<f>) is presented in units of p , using O = and the + sign in (f42|) . 



6 Numerical analysis and the comparison with 
the observational data 

In this section we perform a numerical analysis in order to find constraints on 
the model parameters from the comparison with the observational data. We use 
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the measurements of the Hubble parameter at different redshifts obtained from 
the passively evolving galaxies [53] and baryon acoustic oscillations (BAO). The 
principal idea of the former measurement technique is that for a pair of passively 
evolving galaxies at close redshifts, the differential of their age can be determined. 
Then from the expression H(z) = — Yqrj^f the value of H(z) can be calculated. 
This method has been recently successfully used to constrain the parameters of 
the generalized Chaplygin gas [18] . To constrain the model parameters a and 7, 
we use the H(z) data reported in [54"] . 

We restrict our analysis to the spatially flat Robertson- Walker metric motivated 
by the inflationary expansion and data on CMB anisotropies [55]. Furthermore, 
we fix the share of baryons in the present energy density to = 0.0458 [55], in 
accordance with primordial nucleosynthesis requirements. The expression for the 
Hubble parameter then becomes 

H 2 = S ^f( Pb + p), (44) 

where p b denotes the energy density of baryons evolving with redshift as pb{z) = 
P6,o(l + z ) 3 an d p denotes the energy density of the unified DM-DE component 
evolving with redshift as p(z) = pof(z), where (fTUj) gives 



f(z) 



a + (-wo) 



1/(1-7) 



(45) 



a + (-w(z)) 1 -^' 
The expression for the Hubble parameter then becomes 

H(z) = H (ng(l + zf + (1 - n° b )f(z)) 1/2 , (46) 

where Hq = lOO/i km s _1 Mpc -1 is the present value of the Hubble parameter. 
The functions H(z)/H$ for different values of (ce,7) along with the observational 
data are presented in Fig. HI The x 2 function is calculated according to the 
expression 

2 ^ (H th (z t ) - H° bs (z t ,)) 2 



X 



£ {£1 W-J 1 ™> , (47) 



with the observational data for H(z) taken from [54] . Finally, from (|47|) the prob- 
ability is calculated as 

p = Ae~ x2/2 , (48) 

where the symbol A denotes the normalization. 

Although the theoretically preferred region for a is (0, 1), we run our model 
for a broader range of parameters —5 < a < 5 and 1 < 7 < 10. While we 
marginalize over the parameter h, we select wq = —£1^/(1 — 0,®) = —0.76 (where 
and refer to percentages of the cosmological constant and the baryonic matter 
in the ACDM model [55]) to obtain constraints on the parameters (a, 7) which 
produce the present total EOS consistent with the value of the ACDM model. The 
(68.3%, 95.4%, 99.7%) contours of the marginalized probability density are given in 
Fig EJ They are calculated as curves for which Ax 2 = (2.30, 6.17, 11.8) where Ax 2 
denotes the difference of x 2 a t some point and its minimal value [56} 157] . From 
the figure it is clear that the allowed interval of a grows with the increase of 7. 
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Figure 4: The redshift dependence of the Hubble parameter for several values of (a, 7) 
parameters with = 0.0458 and wq = —0.76. The observed values for the Hubble pa- 
rameters are taken from [54] and the value for the present value of the Hubble parameter 
is H = 70.2 ± 1.4 km s" 1 Mpc -1 [53] . 
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Figure 5: The 68.3% (white), 95.4% (light gray) and 99.7% (gray) probability contours 
in the (a, 7) plane after marginalizing over h. 
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The contour plot for the (68.3%, 95.4%, 99.7%) probability intervals when no 
marginalization over h is performed and the value for the Hubble parameter is 
taken from [55] to be h = 0.702 is presented in Fig. [UJ From figures [5J and E] it is 
evident that the probability contours for the (a, 7) parameters with and without 
marginalization are quite similar. 




-4 -2 2 4 

a 

Figure 6: The 68.3% (white), 95.4% (light gray) and 99.7% (gray) probability contours 
in the (a, 7) plane for h = 0.702 [55] . 



7 Discussion and conclusions 

The analysis presented in the preceding section serves primarily as an orientation 
regarding the allowed part of the parametric (a, 7) space. A more complete analysis 
should consider the growth of inhomogeneities and make the comparison of the 
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model predictions with other observational data, such as the data on the matter 
power spectrum, supernovae of the type la and cosmic microwave background. 
Although preliminary analyses in this direction have already been made [S5], a 
detailed analysis of observational data is left for future work [59J . 

From Figs and it is evident that for larger values of the exponent 7, a 
wider interval of the coefficient a is allowed. A similar conclusion follows from 
the preliminary analysis of the matter power spectrum [S5]. This feature can be 
understood on qualitative grounds using the following argumentation. For large 
values of 7, as long as the EOS parameter is not far from w = 0, the speed of 
sound of the dark component remains small and suppressed by the large exponent 
7. This fact prevents the formation of the sound horizon and its adverse effects on 
structure formation. 

The generalized model of unification of dark matter and dark energy introduced 
in this paper opens a novel perspective on Chaplygin gas and its modifications and 
generalizations. The model defined by the relation (? s = a(— w) 7 encompasses both 
the Chaplygin gas (for 7 = 1, a = 1) and the generalized Chaplygin gas (for 7 = 1) 
as a specific subclass. It could therefore serve as a wider framework for the analysis 
how much these models need to be extended to satisfy the constraint from the ob- 
servational data. From the modeling perspective, the crucial element of our model 
is the relation between the speed of sound, Cg, and the parameter of the EOS, w. 
This relation connects the quantity governing the growth of inhomogeneities with 
the quantity determining the global evolution of the universe. This feature might 
allow easier and more direct transformation of the phenomenological knowledge 
acquired from the data into workable models of the dark sector. In particular, 
an adaptive model assuming piecewise constant values of the speed of sound in 
consecutive redshift intervals is presented in the Appendix. 

A particular challenge for the future research is finding a microscopic explana- 
tion of the dependence c s (w). Here the corresponding microscopic models for the 
Chaplygin gas might serve as a good starting point. In particular, the method of 
representing the evolution in terms of minimally coupled scalar fields, presented in 
section [5] could provide useful information on such microscopic models. 
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useful comments on the manuscript. This work was supported by the Ministry of 
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Appendix 

In the comparison of models where the sound speed of the cosmic fluid is tested 
against the observational data, the goal is to test as general functional form c|(u>) 
as possible. As this is in general quite a formidable task, a powerful approximation 
to the general case is to assume that c| is piecewise constant, i.e. that for every 
interval of redshift [zi-i,Zi), with zq = 0, the speed of sound is constant, i.e. 
Cg = In this approach the parameters are borders of the intervals Zi and 
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the constant sound speed values cr|j. An example of this simple, but powerful 
parametrization was given in [60J for the case of w(z) function. The equations 0, 
(fTU|) and (fTT]) in the redshift interval become 



w(z) 



2 



l+%- 



p(z) = pi_ v 



1 + 



1+I0i_l 



1+z 



3(1+0 



and 



p(z) = Pi_r 



1 + c? 



Wi-l 



1 + 



l+Wi-l / 1 + z \ 3 ( 1+C U 



C Z s i - Wi-l \1 + Zi-l 



2 1 + Wi-l ( 1 + 2 



°s,i 



3(1+^,,) 



(49) 



(50) 



(51) 



The relations between the parameters Wi and pi in neighboring intervals are given 
by relations 

2 / \ 3(1+4,) 



1 



3(1+^,) 



(52) 



and 



Wi-l 



Pi-V 



l + ci 



1 + 



1 + Wi-l ( 1 + z 



Wi-l V 1 + z i-l 



3(1+^,) 



(53) 



The expression (|50|) can be used, e.g. for the comparison of the model with 
piecewise constant speed of sound against the SN la data. 
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